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A self-similar hierarchical solution that is both dynamically and evolutionarily stable is found to 
the multi dimensional Lotka- Volterra equation with a single chain of prey-predator relations. This 
gives a simple and natural explanation to the key features of hierarchical ecosystems, such as its 
ubiquity, pyramidal population distribution, and higher aggressiveness among higher trophic levels. 
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From bacterial colony to human society, hierarchical 
structure is one of the most universal features of ecolog- 
ical systems. It is typically realized as a trophic pyramid 
with exponentially larger populations for lower trophic 
levels. The ubiquity of this hierarchy in nature suggests 
the existence of a simple and robust mechanism behind 
it. 

Let us consider an ecosystem of N species whose pop- 
ulations xi(t), X2(t), Xff(t) are described by a set of 
coupled differential equations with parameters that rep- 
resent the environmental conditions Q, Q ■ Then, the ro- 
bustness of ecological structures is expressed as the pres- 
ence and stability of fixed point solutions. The stabil- 
ity should be considered on two levels of different time- 
scales. The short-term dynamical stability is the stability 
against the perturbation in the dynamical variables x n , 
while the long-term stability is related to the robustness 
of the solutions against the secular variation of environ- 
mental parameters. When some environmental parame- 
ters are at the disposal of a certain species, natural se- 
lection will lead to the realization of the parameter value 
that maximizes the population of that species. This is 
the concept of evolutionary stability jj, |j] . 

There are several numerical simulations of multi- 
species ecosystems that show the spontaneous emergence 
of multi-trophic structure 0, 0. The models employ 
coupled differential equations with stochastic parameter 
variation subjected to evolutionary selection rules. These 
results indicate that pyramidal hierarchy is an evolution- 
arily stable configuration of ecosystems irrespective to 
the fine detail of the model. It is high time to search 
for a simple and clear theoretical explanation before fur- 
ther numerical calculations with increasingly "realistic" 
settings are to be pursued. 

In this article, we consider ecosystems modeled by the 
Lotka- Volterra equation describing N species that form 
a single vertical chain of prey- predator relations. We in- 
tend to prove the existence of hierarchical solutions that 
are stable both dynamically and evolutionarily. 

Let us begin with the two species prey-predator Lotka- 
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Volterra equation 

x\ 

X2 



bxi 



P2X\X 2 



(1) 



-d 2 X 2 + f 2 P2X\X 2 . 



Here b is the reproduction rate of the prey x\ and a the 
environmental limitation factor to its growth. The pa- 
rameter d 2 is the decay rate of the predator x 2 which 
will not subsist without preying on x\ with the rate of 
aggression p 2 . The factor f 2 represents the combination 
of the efficiency of the predation and the average mass 
ratio between prey and predator individuals. All param- 
eters are positive real numbers. Although the dynamical 
variables treated as continuous quantities here, 

they are approximations of the actual integer popula- 
tions. Moreover, in real life, there is a threshold number 
for a population under which a species is not viable. As 
is immediately identified, Q has a nontrivial fixed point 
solution x\(t) = X\, x 2 (t) — X 2 where X\ and X 2 satisfy 



aXi — p 2 X 2 
-d 2 + f 2 p 2 X 1 



(2) 



The eigenvalues A of the linearized map around the fixed 
point are given by 



A 



1 



Ad 2 q(bq - 1 
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(3) 



Therefore, the fixed point is dynamically stable when we 
have f 2 p 2 b > ad 2 . Since, in this work, we are primar- 
ily concerned with the "populous" regime, 6 > a, this 
condition is almost always satisfied. 

Let us now assume that the aggression parameter p 2 
is a quantity that is at the disposal of the predator x 2 
through a long term "behavioral change" . That is, we 
regard Xi and X 2 as functions of p 2 . Naturally, a change 
in p 2 would be directed toward the maximization of the 
predator population x 2 through evolutionary selection. 
We further assume that the time scale for the change 
of p 2 is substantially larger than the time scale for the 
variation of x\(t) and x 2 (t). Then, irrespective to the 
precise mechanism of the variation of p 2 , one eventually 
ends up with the value p 2 that maximizes X 2 (p 2 ). With 
the notation X* = Xi^p^), we have 
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which represents the evolutionarily stable solution. The 
solution is always dynamically stable, since, at these val- 
ues, the real part of the eigenvalue of linearized map 
never becomes positive. The stability against parametric 
variation of p 2 can be judged by 



d 2 X* 
dp 2 2 



P2 



(5) 



Our result shows that the optimal aggression rate from 
the stand point of the predator is to hunt the prey down 
to one half of its natural stability point b/a that is 
reached by x\ when left alone. Already at this point, 
@ gives us some insights. When b and d 2 are compara- 
ble quantities, the predator population X 2 is suppressed 
by the factor / 2 /2 compared to the prey X*. Since / 2 is 
typically smaller than 1, we tend to have a small number 
of predators supported by a large pool of prey biomass as 
a stable configuration. Another interesting point is that 
the milder environment signified by a higher value of b/a 
will increase both X± and X 2 while reducing the opti- 
mal aggression rate p 2 of the predator. This principle of 
noblesse oblige is a widely observed, but nonetheless non- 
trivial aspect of life. This is corroborated, for example, 
by recent field work observation on slave-making ants • 

According to Q), the only way for the prey x\ to in- 
crease its equilibrium population is to "improve the en- 
vironment" by increasing b/a, when there is any such 
mean available to it. An intriguing fact is that "improv- 
ing defensive shield" by decreasing f 2 will not benefit x\ 
directly; it simply decreases the predator population X 2 . 
However, when A 2 is close to the viability threshold, a 
decreasing f 2 would be a sensible strategy for xx, since 
that could drive x 2 out of existence, which would result 
in the instant doubling of 

Next, we consider the case of N = 3 species that forms 
a single chain of prey-predator relations ||| : 



Xx 
%3 



bxx — ax 1 — p 2 XxX 2 , 

-d 2 X 2 + f 2 p 2 XiX 2 - P3X2X3, 
~d 3 X 3 + f 3 P3X2X3- 



(6) 



The fixed point solution is obtained as 

b-aXx- P2X2 = 0, 
-d 2 + hpiXx - P3X3 = 0, 
-d 3 + f 3 p 3 X 2 = 0. 

By rearranging the first two equations, we obtain 



(7) 



&2 — a 2 -^ 2 — P3X3 
-d 3 + /3P3A12 



with 



« 2 



0. 
0. 



do 



(8) 



(9) 



The problem is therefore reduced to the N = 2 case with 
the predator X3 and effectively self-sustaining prey x 2 



which has reproduction and limiting coefficients b 2 and 
a 2 . If the top predator, driven by evolutionary selection, 
tries to maximize its equilibrium population X3 by vary- 
ing P3, it will reach the optimum given by 
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P3 



2a 2 d 3 



x (*) = hh 
3 4a 2 d, 
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*W = X (10) 
2a 2 



The fact that these values are optimum only with a given 
p 2 is indicated by the bracketed asterisk. From the rela- 
tion p^X^ = b 2 /2, one can rewrite the first two equa- 
tions of as 



2 a 



a ( Xx - — I - p 2 X 2 
2 +hP2 [ Xl -Ya 



(11) 



This is essentially the same relationship as in the N = 2 
case ©, with an extra factor 1/2 in front of the first 
terms, and the shift in X\ in the socond. We should 
now suppose that the middle predator x 2 will, in a long 
run, adjust its aggression rate p 2 toward the prey x\ and 
maximize X 2 . We then obtain the solution 



P2 = 



2ad 2 



X* = 



which in turn yields 



8ad 2 d3 



P3 



f 2 b 2 
8ad 2 



x; = 



xi = 3 A 

1 4a' 
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(12) 



(13) 



Note the fact that p 2 here is identical to the N — 2 case. 
We also obtain parametric stability measures as 



d 2 X* 
dpi 



2pf 



d 2 Xj 
dpi 



= —k> ( 14 ) 
P3 



which indeed prove the evolutionary stability of the solu- 
tion. In Fig.l, the phase space profile of one such example 
of evolutionarily stable N = 3 system is depicted. With 
consideration at each stage (jSJ) and §H\ , it is easy to see 
that this evolutionarily stable solution is also dynami- 
cally stable for all parameter values. In effect, the single- 
chain N = 3 Lotka-Volterra equation is broken into two 
N = 2 equations with essentially the same structure, al- 
beit with an additional factor for the lower chain. 

When d 2 and d 3 are comparable quantities, the popu- 
lation of the top trophic level A3 is inherently suppressed 
by the factor compared to that of X 2 , giving a pyra- 
midal profile to the trophic structure. It is amusing to 
note that, from the stand point of the lowest trophic 
species, an N = 3 system, in which two thirds of its 
natural population is left alive, is considerably more "be- 
nign" than an N = 2 system. 

The preceding proof for the N = 3 solution suggests its 
generalization to arbitrary N. This is achieved through 
the realization that the fixed point equation for any 
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FIG. 1: The phase space profile of a N = 3 system at evo- 
lutionarily stable parameter value. The parameters are set 
to be b = 1, a = 1/10000, d 2 = 1/3, d s = 1/5, f 2 = 1/10, 
fz = 1/5. The optimal aggression rates are calculated as 
p2 = 1/1500 and pjj = 1/375. All orbits approach the fixed 
point (Xt,X%,X£) = (7500,375,64.5). 



mid-level population x n can have both purely prey-like 
and purely predator-like representations. Let us start 
with the N vertically-coupled Lotka-Volterra equation 
with evolutionarily adjustable aggression parameter p n 
for each species x n 



(15) 



X\ — bxi — ax\ — P2X1X2, 

X 2 = -d 2 X2 + j2P2X\X 2 - P3X2X3, 

%n — d n X n -\- f n p n X n — \X n p n ^-\X n X n -]-i , 



x N = -d N x N + f N pNX N -iX N . 

In general, a trophic level can comprise several compet- 
ing species. In our simplified treatment, however, such 
species are lumped into a single population variable. The 
equations for the nontrivial fixed point (Xi, Xn) are 

b-aX 1 -p 2 X 2 = 0, (16) 
-d 2 + /2P2X1 - p 3 X 3 = 0, 



~d n + /nPn^n-l — Pn+\X n +\ — 0, 



— dN + JnPnXn-i — 0. 

Apart from the species with the highest trophic level Xn, 
each of these can be transformed to the form 



b n — a n X n — p n+ iX n+ i — 

with the recursive definition 
1 



p 2. 

a n — JnPn" 



a n -i 



f bn-1 _ , 
0>n— 1 



(17) 
(18) 



Rewriting (|17fl . we obtain the "slave" form 



r\nb n - a n [ X n — (1 - r] n ) — ) - p„+iA"„ + i = (19) 



for n = 1, .., N — 1. Let us now assume the relation 

(*) . VnK 

^n+l n+l — 2 ' 

If we combine l|16fl. I|20|) with a requirement 

Vn+l 



(20) 



1 - Vn 



we obtain the "master" form 



(21) 



- Vndn+l + fn+lPn+1 X n - (1 - T) n ) = (22) 

\ a n J 

for n = 1 , . . , N— 1 . The equations l|16l) are now decoupled 
to (JV — 1) pairs of prey-predator equations i|19|) and l|22|l . 
We then have 



and 



(*) _ 2a n rf„ + i 
Pn+i — — f 7 ! 

Jn+lOn 



.(*) _ fn+lbl „^ , b 



(23) 



X nil=Vnf^, xW = (2- % )-^. (24) 

This result justifies the assumption (|20l) a posteriori, and 
the whole procedure becomes consistent. From the last 
equation of (|16(l . we observe that ?7Ar_i should be set to 
one, which results in i]n-2 = 1/2, = 3/4, • ■ ■ . We 

finally obtain the following explicit forms for the evolu- 
tionarily and dynamically stable solution: 



x! = 
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B N b 
2 N ~ 1 a 

BN-lf2b 2 

2 N ad 2 ' 



(25) 



P2 



2ao?2 



B 1 f 2 ---f N b 2 



2™-*ad 3 ' Pn f 2 ---f N b 2 
The stability with respect to the variation p n is given by 

B]^- 2 b 



d 2 Xj 
dp 2 2 



Bjsi-ib d 2 X£ 



dp 2 



d 2 X* N 



3' 



(26) 



d P N P N 

Here the coefficient B n is a variant of the Fibonacci series 
defined by 



B 



n+2 



B n +i + %B n , Bi — B 2 — 1. 



(27) 



Some of the numbers are B 3 — 3, B4 = 5, B5 = 11, 
B 6 = 21, B 7 = 43, • • • . 
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TABLE I: The evolutionarily stable hierarchical population for N species Lotka-Volterra equation up to N 
population of n-th trophic level, and p n its aggression rate toward its prey. 



5. X* is the 
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N = 4 



Pn 



.f2/ 3 / 4 6 2 

64ao?4 
/ 2 / 3 ° 2 



32ac?3(j4 

8ac?2^ 3 
/ 2 / 3 fo 2 

2ad2 
hb 



1 D- 



56 

8a 



AT = 5 
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3/ 2 / 3 & 2 
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In table I, the hierarchical solutions up to = 5 are 
listed. The most notable feature is of course the expo- 
nentially smaller population in higher trophic levels. As- 
suming d n ~ b for all n, we have a decrease in the pop- 
ulation by factor f n /2 for each increase of one trophic 
level. Since f n is in general substantially smaller than 
one, we get a pyramidal hierarchy with a steep exponen- 
tial decrease. We should also mention the self-similarity 
of the solution: For any given trophic level, the portion 
of its "natural" population saved from exploitation by 
higher trophic levels varies like 1/2, 3/4, 5/8, • • • , when- 
ever more trophic levels are added on top. On the other 
hand, its optimal aggression rate is unaffected by the 
presence of higher trophic levels. A higher value of p* n 
for larger n is a direct result of the scarcity of its prey. 

Ultimately, the quantity b/a gives the base biomass, 
on top of which the whole trophic pyramid structure is 



built. Since there is a minimum population for the high- 
est trophic species to be viable, this naturally puts a limit 
to the maximum number for the trophic hierarchy of an 
ecosystem with a given base biomass. 

In summary, within the framework of a single vertical 
food chain model, a pyramidal self-similar hierarchy is 
found in Lotka-Volterra system. It might be possible 
to generalize our results to models with plural species 
in each trophic levels. Hopefully, the search for generic 
properties of Lotka-Volterra system along this line shall 
provide a solid backbone for experimental and numerical 
studies of ecosystems. 
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